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Particle Acceleration Very Near an X-Line in a Collisionless Plasma 


L. R. LYONS AND D. C. PRIDMORE-BROWN 

Space and Environment Technology Center , The Aerospace Corporation, 
El Segundo, California 90245 


In a previous paper, we applied a simplified model for particle motion in the vicinity of a mag- 
netic X-line that had been introduced by Dungey. We used the model to quantitatively show 
that an electric force along an X-line can be balanced by the gyroviscous force associated with 
the off-diagonal elements of the pressure tensor. Distribution functions near the X-line were 
shown to be skewed in azimuth about the magnetic field and to include panicles accelerated to 
very high energies. In the present paper, we apply the previous model and use the distribution 
functions to evaluate the energization that results from particle interactions with the X-line. 
We find that, in general, this interaction gives a spectrum of energized particles that can be 
represented by a Maxwellian distribution. A power-law, high-energy tail does not develop. 
The thermal energy, K, of the Maxwellian can be expressed simply in terms of the field parame- 
ters and particle mass and charge. It is independent of the thermal energy, of the particle 
distribution incident upon the region of the X-line, provided that Kj < K. Significant energiza- 


tion is not found for Kj > K. 

INTRODUCTION 

Energy transfer in space plasmas is often associated with 
a magnetic X-line along which lies an electric field [e.g. 
Vasyliunas, 1975], and the physics of the forces and particle 
acceleration near such X-lines is an active research topic. In 
this paper, we consider the acceleration of particles very near 
an X-line using particle distributions that are consistent with 
the requirement for force balance near the X-line 

Plasma dynamics near an X-line are often addressed 
within the framework of MHD under the assumption that 
some sort of resistivity maintains the force that opposes the 
electric force along the X-line. However, space plasmas are 
generally collisionless so that the applicability of the con- 
cept of resistivity is highly questionable. 

It has been pointed out [Vasyliunas, 1975; Sonnerup, 
1979, 1988] that force balance along an X-line in a colli- 
sionless plasma requires a balance between the electric force 
and the gradient of the off-diagonal elements of the electron 
pressure tensor (referred to as "gyroviscosity"); however, 
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identification of a process that could lead to such gyroviscos- 
ity was not accomplished until recently. The breakthrough 
on this topic was by Dungey [1988], who evaluated colli- 
sionless particle motion near an X-line using a simplified 
model. He showed that the particle motion near an X-line 
will lead to a skewing of the velocity-space particle distribu- 
tions and thus to significant off-diagonal pressure tensor 
elements. He suggested that the resulting gyroviscosity 
could be sufficient to balance the electric force along the X- 
line. Such gyroviscosity should be present in simulations 
of collisionless reconnection and should be identifiable by 
evaluation of the off-diagonal elements of the pressure 
tensor. 

We [Lyons and Pridmore-Brown, 1990; hereafter referred 
to as paper 1] applied Dungey's model to evaluate particle 
distribution functions very near an X-line, and we used these 
distributions to quantitatively evaluate the force balance. 
We found that, within the limits imposed by Dungey's 
model, gyroviscosity can balance the electric force along an 
X-line in a collisionless plasma independent of other plasma 
and field parameters. The distribution functions we evalu- 
ated include the effects of particle energization from electric 
field acceleration in the weak magnetic field region near the 
X-line. Here we use these distributions to evaluate this 
energization. 
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DUNGEY* S MODEL 

The magnetic field (B) model used to evaluate the particle 
trajectories and distribution functions is two-dimensional, 
having d/d y = By = 0. The X-line is taken to lie along the 
y axis and to be imbedded within a current sheet that is cen- 
tered in the z = 0 plane. B x reverses in sign across the cur- 
rent sheet Symmetry is assumed about the X-line, so that 
B x (z) = -B x (-z), and B z (x) = -B z (-x). The normal component 
of B across the current sheet is expanded in a Taylor series 
about the X-line, so that B z = Px. The electric field Ey is 
taken to be uniform and to be directed parallel to the current 
in the positive y-direction. 

The basic motion of particles in the vicinity of a current 
sheet having a weak but uniform B z was originally presented 
by Speiser [1965] and is described in paper 1. Particles that 
reach the midplane of the current sheet oscillate about the 
midplane and simultaneously gyrate about B z for one-half a 
gyro-circle; the gyration begins with an x-component of 
velocity, v X o» directed towards the X-line, and v y0 = 0. 
This gives a current in the positive y-direction, as is required 
to maintain the current of the current sheet. Only particles 
with qv y > 0 remain in the current sheet. This is because 
particles with qv y < 0 are accelerated away from the current 
sheet by the magnetic force due to B x (which the electrons 
encounter as a result of their oscillatory motion about z = 
0 ). 

Following Dungey [1988] and paper 1, we take advantage 
of the decoupling that occurs between the x-y component 
and the z component of motion for particles as they oscillate 
about the z = 0 plane [Speiser, 1965, 1968]. This decou- 
pling is assumed to hold as the x-y component of a particle 
trajectory crosses an X-line. Speiser did not specifically 
include this situation or the effects of 3B z /dx * 0 in his 
analysis, though particle behavior near an X-line has been 
considered by Stem [1979]. The evaluation of particle 
motion within the midplane of the current sheet can thus be 
simplified by calculating only the x-y component of the par- 
ticle trajectory at z = 0. This motion is taken to begin with 

Vy = 0. 

Under the above conditions, the equations of motion for 
particles of charge q and mass m are: 

m(dv x /dt) = qpxvy (1) 

m(dvy/dt) = q(-pxv x + E y ). (2) 

Using the normalized temporal and spatial variables intro- 
duced by Dungey [1988], 

x' = x/L, where L = IE y m / (p 2 q)l 1/3 (3) 

i = t/t, where t = L 2 p / IE y l , (4) 


the equations of motion become 

dv x 7dt’ = x'vy'sign(q) (5) 

dvy'/dt’ = (-x’v x ‘ ± l)sign(q), (6) 

where ± refers to the sign of Ey, which we take to be posi- 
tive in the present work. Note that Ey, m, (i, and q all have 
a magnitude of unity in normalized units. Also, since the 
equations of motion are applicable to particles of any charge, 
we take sign(q) = +1. Thus, our results are applicable to 
electrons and ions provided the sign of vy' is changed for 
electrons. 

Representative particle trajectories in the vicinity of the 
X-line at x' = 0, obtained from solving equations (5) and (6) 
using a Runge-Kutta scheme, are shown in Figure 1. The 
trajectories shown are for particles that start their motion 
within the current sheet at values of x' > 0. The distortion 
of the trajectories from semicircles results from the accelera- 
tion by the electric field and the variation of B z with x\ 

VELOCITY SPACE DISTRIBUTIONS 

Trajectories such as shown in Figure 1 can be used to 
evaluate normalized distribution functions, f , in phase space 
for the z’ = 0 plane if a distribution function, f|', is specified 
for particles incident upon that plane. In paper 1, we illus- 
trated the features of these distributions by taking a "top hat” 
form for fi*. The value of f|’ was taken to be a constant for 
all initial velocities, v x0 \ directed toward the X-line and 
having 0 < lv xo ’I ^ 1, and fi' was taken to be zero otherwise. 
Particles were assumed to be incident at all values of x\ 
With this incident distribution, the distribution at all x 1 in 
the z’ = 0 plane can be described by an area in v x ,v y ’ space 
within which f is constant and outside of which f = 0. It 
thus suffices to trace the boundary of this area in order to 
describe f (v x \v y *) for an x' at z' = 0. 

Boundaries in velocity space within which f is constant, 
as obtained in paper 1 for selected values of x\ are shown in 
Figure 2. The boundaries show that f becomes increasingly 
skewed towards positive v x and v y ' as x' 1. This is the 
skewing that gives rise to positive values for the off- 
diagonal elements of the pressure tensor and thus to the 
gyroviscosity. In paper 1, we used these distributions to 
evaluate moments of f and the balance of forces near the X- 
line. 

While all particles used to construct the distributions in 
Figure 2 entered the z' = 0 plane with Iv'l < 1, the distribu- 
tions show tails extending to very high values of v\ These 
tails result from particles having trajectories that become 
nearly parallel to the X-line at very small values of x' (see 
trajectories in Figure 4 of paper 1), and the tails are associ- 
ated with a net energization of particles that results from 
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Fig. 1. Representative particle trajectories in the x’-y 1 plane 
near an X-line at x' = 0. The sign of y' is correct for positively 
charged particles. For negatively charged particles, the sign of 
y* must be changed. 

their interaction with the X-line. 

To evaluate the energization of particles from their inter- 
action with the X-line, we assume a Maxwellian, 

Fi'(vxo') = exp[-v X o' 2 /vi’ 2 ] , (7) 

for the distribution of initial speeds of particles following 
their incidence upon the z = 0 plane. Under the assumption 
that initial velocities are directed towards the X-line, equa- 
tion (7) corresponds to an incident distribution function in 
velocity space given by: 

fiXvxo’.Vyo’.vzo^ = exp[- Vxo^/Vi* 2 ] Stvyo'jfiXvzo'), (8) 

for xVxo' ^ 0, 

where JffCvzoldvzo* = 1, and fi' = 0 for xv X o’ > 0* 
Applying the assumption that the z and x-y components of 
particle motion are decoupled within the z = 0 plane, we take 
(vz^dvz' = 1 throughout the z’ = 0 plane and ignore the v z ’ 
dependence of f . 

Using equations (3) and (4), the normalized thermal 
velocity of incident particles can be expressed in terms of 
real parameters as 

Vi’ = [2(Kiq 1/3 P 2y3 )/(Ey 4/3 m 1/3 )] 1/2 , (9) 


where K| is the thermal energy of the distribution expressed 
in eV. If we chose as nominal parameters for the distant 
geomagnetic tail Kj = 110 eV, P = 0.3 nT/Re, and E y = 
0.25 mV/m, then equation (9) gives vf = 1.0 and 0.286 for 
electrons and protons, respectively. 

With a Maxwellian distribution for FjXvxo’), f(vx\vy') 
does not divide into regions where f = constant and where f = 
0. However, the regions of maximum f remain similar in 
shape to the boundaries shown in Figure 2, so that the dis- 
tributions in Figure 2 still qualitatively represent frvx’.vy 1 ). 
To quantitatively evaluate f(v x \vy') at a given x\ we use 
equations (5) and (6) to trace particle trajectories backwards 
in time from points in v x \vy* space until v y ' = 0. In this 
way, we find initial speeds, v X o> corresponding to the 
points in vx\vy* space, and we obtain f from equation (8). 

Since our goal is to evaluate the energization from parti- 
cle interactions with an X-line, we choose values of v' = 
( v x 2 + Vy' 2 )^ and determine f over a range of values for f 
= tan“l(vy7v x r ) selected so as to allow the evaluation of 
F(v f ) = Jf(v»d<|>/n. This gives the distribution of speeds 
for particles within the current sheet after their interaction 
with the X-line, which can be compared direcdy with the 
incident distribution given by equation (7). As can be seen 
from Figure 2, the range of important $ values becomes 
increasingly narrow with increasing v\ and, in fact, becomes 
several orders of magnitude less than 1° wide for the highest 
v' values of interest. This narrow range would be difficult to 
handle in a model that includes three-dimensions of particle 
motion; however, it does not present any problems here. 

Comparisons between F and Ff, plotted versus v'/vf for 
x 1 = 1, are displayed in Figure 3. In principal, the results are 
only dependent on vf . However, we show results for values 
of vf that correspond to specific real parameters. Results for 
the values of vf given above that correspond to electrons and 
ions for nominal tail parameters are shown in the middle 
panel (identified as Ey = 0.25 mV/m). Results are also 
shown for values of vf that have been adjusted, using equa- 
tion (9), for Ey « 0.062 mV/m and for Ey = 1.0 mV/m 
with Ki and p unchanged. In each panel, the solid line gives 
Ff. When plotted versus v’/vf, Ff is the same for all vi', 
which facilitates comparison of the amount of energization 
between the different cases. Our evaluation of F(v’) for each 
Ey is given by filled circles for the value of vf correspond- 
ing to electrons and by filled triangles for the value of vf 
corresponding to protons. 

Figure 3 shows significant energization for all cases hav- 
ing vf < 1, the energization increasing with decreasing vf. 
Such energization does not occur for vf 2 » 1, as can be seen 
for electrons with Ey = 0.062 mV/m, which corresponds to 
vf = 2.5. Keeping other parameters fixed, the energization 
is greater for protons than for electrons, since vf « nr 1/6, 
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Fig. 2. Boundaries in normalized velocity space within which f is constant at selected values of x' as obtained in paper 1 for 
the "top-hat" distribution of particles incident upon the current sheet. 


and it increases with Ey since vj’ «= Ey'^. 

The distributions in Figure 2 suggest that the energy dis- 
tributions few* the particles might show a high energy tail 
since some low-energy particles are accelerated to very high 
energies. However, Figure 3 shows that F(v') falls off as a 
Maxwellian for v 1 > vi’. In fact, we have found that F(v’) at 
v' > 1 is very well represented by a Maxwellian having a 
normalized thermal velocity equal to two for all the cases 
examined that have < 1. Specifically, the Maxwellian 

F(v') = 0.5 v’3/2exp[-v72]2 ( 10 ) 

which is shown in Figure 3 by dashed lines, can be seen to 
fit our results very well for v* > 1 . All our distributions 
have F(v*) — > 0.5 as v' 0; thus equation (10) does not fit 
F(v') for v' « 1. 

The results in Figure 3 are all for x 1 = 1. Results for the 
nominal tail parameters are shown in Figure 4 for x' = 0, 1, 
and 2, which show that F(v r ) is essentially independent of x' 
near the X-line. 


The results in Figures 3 and 4 show that the energy dis- 
tribution of particles in a current sheet very near an X-line 
should have a thermal velocity of approximately two in 
normalized units, independent of the thermal energy of the 
particles incident upon the current sheet, provided the nor- 
malized thermal velocity of the incident particles vi’ < 1 . 
This implies, from equations (3) and (4), that such particles 
will have a thermal energy K in eV given by 

K = 2E y 4 /3ml/3/ql/3p2/3 (11) 

independent of Kj provided K > 4Kj. No energization 
occurs if equation (11) gives K < Kj. For our nominal tail 
parameters (Ey = OJ25 mV/m), equation (1 1) give a thermal 
energy of 0.4 keV for electrons and 5.3 keV for protons. 
These are of the order of the thermal energies of particles 
observed within the tail plasma sheet; however, the effects 
of particle motion away from the X-line region would need 
to be evaluated for a direct comparison with plasma sheet 
observations to be meaningful. 
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V / Vj’ 



Fig. 3. Normalized distribution functions F’(v’) versus v'/Vj'. 
In each panel, our evaluations of F at x' = 1 are given by filled 
circles for electrons and by filled triangles for protons. The 
Maxwellian approximation for F(v’), equation (10), is given 
by the dashed lines, and the distribution of incident particles, 
equation (7), is given by the solid lines. The values of v^ used 
for electrons and ions, respectively, are 2.50 and 0.715 in the 
upper panel, 1.00 and 0.286 in the middle panel, and 0.400 and 
0.113 in the bottom panel. These values correspond to the 
electric fields given in each panel and p = 0.3 nT/R^ 

SPECTRA OF PARTICLES EJECTED FROM 
CURRENT SHEET 

The energy spectra calculated in the previous section 
would only be directly observable near the mid-plane of the 
tail current sheet and within a normalized distance x* < 2 
from an X-line. This distance is quite small for electrons. 


Fig. 4. Same as Figure 3, except that results are shown at x' = 
0, 1, and 2 for the values of vj' corresponding to the 0.25 mV/m 
electric field. 

and is given by 2L = 1700 km for our nominal tail parame- 
ters. For protons, 2L = 20,800 km. However, after travers- 
ing half a gyro-orbit about B z , all particles in our model are 
presumed to be ejected from the current sheet These parti- 
cles will flow along field lines just within the boundary 
between open and closed field lines, where they are more 
likely to be observed than within the current sheet. It thus 
is worthwhile to consider the energy spectra of the ejected 
particles. 

Within the context of our model, it is not sufficient to 
simply evaluate the distribution function of ejected particles 
as a function of energy at a specified x\ This is because of 
the strong variation of the location xf of particle ejection as 
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a function of initial velocity v x0 ’ and position Xq\ In par- 
ticular, particles that attain high final velocities Vf (i.e., Vf 

> 3) are ejected over very narrow ranges of xf* that become 
increasingly narrow, but occur at increasing x\ with increas- 
ing vf\ For example, most particles with vf' = 3 are ejected 
within a range -0.2 wide in x' centered at xf* = 2.6, whereas 
particles with vf' = 6 are primarily ejected within a range - 
0.0002 wide in x' centered at xf' = 3.536. (These x' ranges 
have been determined using fj' for v^ =1.) Asa result, an 
ejection spectrum at a specified x* would consist mosdy of 
particles having a very limited range of energies. 

We thus have decided to evaluate the total number of 
ejected particles that have interacted with the region near the 
X-line as a function of vf\ We do this by assuming that 
particles are incident upon the current sheet within a range of 
x 0 * centered about x' = 0, and evaluating FjXvf*) = 
Jf (v^xf'Jdxf 1 for various values of vf\ This integration 
gives us an energy spectra that can be interpreted as the spec- 
tra averaged over the spatial extent of the particles' gyromo- 
tion near the X-line. Here we take the range of incident par- 
ticles to be Ixq*! < 2.25. This particular range was chosen 
because 2.25 is the minimum Ixfl for incident particles hav- 
ing a speed of zero, the speed at which F(v*) maximizes, and 
we thus obtain a continuous spatial distribution for these 
particles. The precise size of this range, however, is not 
important to our results. 

Energization of particles to high final velocities maxi- 
mizes very close to x‘ = 0, whereas these energized particles 
are ejected beyond 1x1 = 2. Thus, at higher speeds (i.e. vf' 

> 2), we expect the distribution of ejected particles FxXvfO 
to be similar to the distributions F(v T ) within the current 
sheet at x' < 2 that are shown in Figures 3 and 4. Figure 5 
shows this similarity for the values of v^ that correspond to 
our nominal tail parameters. In Figure 5, Fx , (vf') is given 
by open circles for electrons and by open triangles for pro- 
tons. The similarity between FxXvf') and F(v’) is also 
found for the other cases shown in Figure 3. (All particles 
are ejected with vf* > 1, so that FjXvf') goes to zero at low 
values of vf*.) Thus, the energy spectra given by equation 
(10) is appropriate both within the current sheet and for the 
particles ejected from the X-line region. Again, however, an 
actual comparison with observations will require additional 
considerations of particle motion. In particular, the effects 
of particle mirroring and the resulting multiple interactions 
with the current sheet need to be evaluated 

COMPARISON WITH PREVIOUS RESULTS 

Several studies have been performed that include the full 
three-dimensions of particle motion near an X-line. In this 
section, we compare our results to results from these stud- 
ies. The comparisons provide a test of the assumptions 



Fig. 5. Same as the middle panel of Figure 3, except the calcu- 
lated distribution FxX v f ) of particles ejected from the current 
sheet has been added. Values of the ejected distribution are 
given by open circles for electrons and by open triangles for 
protons. 

applied in our model to limit the evaluation of particle 
motion to two-dimensions. When three dimensions of 
motion are included, it becomes difficult to evaluate the tra- 
jectories of a sufficient number of particles to determine the 
spectra of energized particles to high energies, where the 
value of F(v) is many orders of magnitudes below its peak 
near v = 0. To the best of our knowledge, such a study has 
not been carried out. Also, none of the three-dimensional 
studies have considered the possibility of significant 
azimuthal asymmetries about B, which lead to the gyrovis- 
cosity and are crucial to the physics near the X-line. 
Nevertheless, a number of interesting results have been 
obtained that are relevant to the present study. 

Martin and Speiser [1988] evaluated particle distributions 
that have been averaged over the azimuthal angle around B. 
They found a ridge of enhanced values of f in velocity space 
that results from particle energization near the X-line in their 
model. This ridge is probably the same feature as the tail 
extending to high energies in Figure 2; however, the skew- 
ing that gives rise to the gyroviscosity cannot be determined 
from their azimuthally averaged distributions. 

Burkhart et al. [1990] considered cold particle distribu- 
tions incident upon the region near an X-line. They found 
that the velocity of the particle distribution flowing awav 
from the X-line was proportional to Ey^ql^/ml'^pl/^ 
which is the same scaling relation that we have obtained 
[equation (11)] for the thermal velocity. They found that the 
outflowing particles were magnetic field-aligned. This 
information cannot be obtained from our evaluation of parti- 
cle motion. However, if we assume that particles ejected 
from the current sheet in our model are field aligned, then we 
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obtain the same scaling as did Burkhart et aL for the outflow 
velocity. Thus, their results provide an excellent verifica- 
tion of the scaling that we have obtained. 

Deeg et al. [1991] followed the three-dimensional trajec- 
tories of 14,000 cold (~1 eV thermal energy) protons near an 
X-line and were able to obtain a statistically significant 
energy spectrum of energized particles for one set of condi- 
tions. They assumed that particles are incident upon the X- 
line over a limited range of x near the X-line, and they 
counted the total number of particles ejected from the X-line 
region as a function of energy. This procedure is essentially 
the same as ours, so that their spectrum provides an excel- 
lent test for our two-dimensional results. Their spectrum is 
shown in Figure 6 as the total number of particles per unit 
energy N(K)dK [which is proportional to v^Fx(v)] versus 
particle energy. 

They used Ey = 0.3 mV/m and p = 2.5 x 10'^ T/m. 
From equation (1 1), our analysis predicts a Maxwellian dis- 
tribution of energized particles having a thermal energy of 
480 eV. The dashed line in Figure 6 shows a Maxwellian 
distribution having a thermal energy of 300 eV (which pro- 
vided a better fit then did 480 eV), which can be seen to fit 
the simulated data very well. This comparison shows very 
good agreement between the general results from our two- 
dimensional model and the three-dimensional results for the 
particular set of parameters used by Deeg et al. [1991], the 
difference in thermal energies being not very significant. 
Also, using 500 particles, Deeg et al. investigated how the 
energization varies with Ey and p, and they found the same 
variation that we have found. 

Note that even with 14,000 particles, Deeg et al.'s [1991] 
spectrum does not extend to high enough energies to distin- 
guish a power-law high-energy tail from a Maxwellian. Our 
spectra do show this distinction, whereas their spectrum at 
energies above - 0.5 keV can be fitted as well by a power 
law as by a Maxwellian. However, the excellent agreement 
we find with their results gives strong support to the valid- 
ity of our general two-dimensional result for the particle 
energization. 

Goldstein et al. [1986] and Ambrosiano et al. [1988] 
evaluated the acceleration near X-lines and included the 
effects of magnetic turbulence obtained from MHD simula- 
tions. They found significant acceleration due to particle 
trapping in magnetic bubbles that developed in the simula- 
tions. Our results, however, are not directly comparable to 
their results because we did not allow for turbulence and 
because the basic physics of collisionless particle motion 
near an X-line (e.g., the gyro viscous force and the associated 
asymmetries in particle distribution functions about B) does 
not exist in MIS). 

Finally, it is interesting to compare the results for parti- 



Fig. 6. Total number of particles per unit energy ejected from 
the region near an X-line as calculated by Deeg et al. [1991] 
from a three-dimensional simulation of particle motion. The 
dashed line gives a Maxwellian distribution for 0.3 keV ther- 
mal energy. 

cle energization near an X-line with those obtained by 
Speiser and Lyons [1984] for a weak, but constant, B z 
across a current sheet. They found that initially cold parti- 
cles were accelerated to a final velocity vf = 2Ey/B z . Our 
present results give a thermal velocity vtb = 
2E y 2/3 q 1/ 3/ m l/3pl?3 = 2E y /B z (x' = 1). Thus the present 
results give a thermal velocity that is the same as the final 
velocity in the B z constant model, if B z is evaluated one 
scale length L away from the X-line. 

CONCLUSIONS 

We have applied Dungey's [1988] two-dimensional model 
of fields and collisionless particle motion to evaluate the 
particle energization near a magnetic X-line. We find that 
this process results in an energy spectrum of particles that 
can be approximated as a Maxwellian. The same 
Maxwellian applies for both the particles within the current 
sheet in the vicinity of the X-line and for the particles 
ejected from that region of the current sheet. This 
Maxwellian has a thermal energy, K, given by equation 
(11), that depends upon field parameters but is independent 
of the thermal energy, Ki, of the particles incident upon the 




170 PARTICLE ACCELERATION NEAR AN X-LINE 


region of the X-line. It is necessary that equation (11) give 
K > Ki for interactions with an X-line to cause energiza- 
tion. No energization occurs if equation (11) gives 
K < K|, and (1 1) is not applicable to such a situation. 

For typical magnetotail parameters, the thermal energy 
given by equation (1 1) is a few hundred eV for electrons and 
a few keV for protons. These energies are of the order of 
magnitude of those observed in the distant tail plasma sheet, 
which suggests that interactions with the neutral line in the 
distant tail may contribute to the energization of plasma 
sheet particles. During geomagnetic disturbances, induction 
by time-varying magnetic fields may yield electric fields in 
the tail that are significantly enhanced over typical values. 
Taking E y = 4 mV/m as perhaps an extreme example, we 
obtain thermal energies of 17 keV for electrons and 210 keV 
for protons. This is consistent with the suggestion that 
interactions with a neutral line may contribute to energetic 
particle phenomena that have been observed in the tail dur- 
ing substorms [e.g., Armstrong and Krimigis , 1968; 
Krimigis and Sarris, 1979; Baker at al., 1979]. 

We have compared our results to results from three- 
dimensional studies of particle motion near an X-line. 
These comparisons show our results to be consistent with 
the parameter scalings and the one detailed energy spectra 
obtained from the three-dimensional studies. Our results 
compliment the three-dimensional studies by giving an accu- 
rate energy spectrum to high energies that can be scaled to 
different magnetic and electric field parameters, and by pro- 
viding details of the particle distribution functions such as 
the azimuthal asymmetries very near an X-line. 
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